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Abstract 



A lattice approach is developed to measure the sphaleron free energy. Its feasibility is 



demonstrated through a Monte Carlo study of the two-dimensional 0(3) sigma model. 



In the electroweak standard model, baryon number (or, more precisely, B + L) is violated 
C*~> ■ by an anomaly |ffl]. Baryon number violation is associated with transitions between different, 

topologically distinct vacua. These vacua are separated by an energy barrier. Sphalerons Q 
are classical solutions corresponding to saddle points on top of the barrier. The lowest-energy 
sphaleron has an energy of the order of 10 TeV (somewhat depending on the Higgs mass). This 
Q-i, energy is concentrated in a region of size rri^. 

(— ! \ 

The tunneling rate between nonequivalent vacua is tiny, due to the small electroweak gauge 
coupling. Therefore, baryon number violation in the standard model was not considered to be 
^ ■ of practical importance for some time. Later, however, it was realized that the energy barrier 

could be overcome more easily by classical transitions instead of tunneling. Classical transitions 
become relevant at high temperature 0], or perhaps in high energy scattering ||, |5|. 

The case of high energy scattering is still badly understood (see e.g. Ref. || for a recent re- 
view). The case of high temperature is extremely interesting because anomalous baryon number 
violation could provide a scenario for the generation of the baryon number of the universe. 

In the high temperature phase with restored symmetry (B+L) transitions are believed to be 
frequent, with a rate [0, ||, ||, 10] 

T = 1 {a w Tf (1) 

per unit volume, where 7 is some constant of order one. This would have the striking consequence 
that any (B+L) asymmetry generated at a scale close to the Planck mass would be washed out 
completely. Where would the baryon asymmetry of the universe come from in that case? 
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One very attractive possibility would be the generation of a (B+L) asymmetry at the elec- 
troweak phase transition [11, 12, [l3|, 14, [H], [l6|| . This would be possible if the electroweak phase 



transition would be strongly first order. Presumably, this requires a more complicated Higgs 
sector (not just one doublet), which would be an interesting prediction. The electroweak phase 
transition is intensively studied by perturbative |l7], [D| as well as lattice techniques ]l9|, ^0| . 

The survival of such a (B+L) asymmetry after the electroweak phase transition is a non- 
trivial matter, however. It could be washed out afterwards by transitions across the sphaleron. 
This will not be the case if the sphaleron barrier near the critical temperature is high enough. 
The corresponding condition on the sphaleron free energy F sp (T) is found by comparing the 
transition rate with the expansion rate of the universe. 

^1 > 45, (2) 

where T c is, strictly speaking, not T cr n, but the temperature at which the phase transition is 
completed. In the case of a second order or weakly first order phase transition this relation will 
not be satisfied. 

The trouble is that the temperature dependence of the sphaleron free energy is not very 
well known, in particular close to the critical temperature. Perturbation theory is plagued 
with severe infrared divergences in this regime. This should also apply to the treatment of 
fluctuations around the sphaleron. The usual way of estimating the sphaleron free energy is by 
rescaling, i.e. replacing the zero-temperature Higgs expectation value by v(T), obtained from 
the temperature dependent effective potential. This is correct to leading order but could be 



quite misleading because of infrared singularities of higher loop contributions [21]. The role of 
the effective potential itself is also obscure because the sphaleron receives its energy mostly from 
the non-convex region which is unphysical. 

Therefore a direct determination of the sphaleron free energy, i.e. the transition rate, is an 
important task. The sphaleron free energy is directly related to the effective potential of the 
Chern-Simons number. A fraction of configurations with CS number close to ^ will contribute to 
the anomalous processes in question. The basic objective is therefore to measure the probability 
of configurations close to the sphaleron, and derive the sphaleron free energy. 

Determining the probability distribution of the CS number requires a reliable way of measur- 
ing it. In a typical Monte Carlo configuration, quantum fluctuations (completely unrelated to 
topological features) may easily add up to produce a CS number of, e.g., 0.5. A straightforward 
measurement of the CS distribution would therefore overestimate the rate of fermion number 
violation. By focussing on the known properties of the sphaleron, it is possible to select con- 
figurations which look like an underlying classical sphaleron with added quantum fluctuations. 
Even so it is necessary to check the stability of such assignments. This situation is well known 
from lattice estimates of the topological susceptibility of pure Yang-Mills theory. In that case, 
cooling techniques p^, S3| have been shown to be useful. Instantons show up as long-lived states 



2 



in this way. This method cannot be taken over without modification to measure the CS density, 
however, because sphalerons are not stable classical solutions like instantons, but saddle points. 
They cannot be found by simply minimizing the action. A different smoothing procedure is 
needed. We propose to use the square of the tadpole (i.e. the square of the equation of motion) 
as a new 'action' density of the cooling procedure. Any classical solution, stable or not, will now 



appear as an attraction point 24]. 



In order to test such an algorithm one should use a simple model, where the magnitude of 
topological transitions is sufficiently understood. The Abelian Higgs model in 1+1 dimensions 
would be a possible candidate. We prefer the study of the 1+1 dimensional 0(3) sigma model 
with some external magnetic field, because of its similarities with the 4-dimensional standard 
model. The action is _ 

-dfj.fi- dfjfi + L; 2 (l + n 3 ) 



S = I d z x 

r 



2 



(3) 



with the constraint 



ft 2 {x) = 1 (4) 

Without the external field uj, the theory is asymptotically free and possesses instanton solu- 
tions. Symmetry breaking due to the external field removes the instantons as true solutions, 
but sphalerons as unstable saddle point solutions appear, just like in the SU(2) Higgs theory. 
Topological transitions in this model have been studied in quite some detail by Mottola and 
Wipf |^| . The one-loop expression for the transition rate has been worked out and is expected 
to be reliable at weak coupling (this is not obvious in the standard model because of infrared 
problems near the phase transition). To exponential accuracy the transition rate is given by 

F S p(T) 

r ~ e t (5) 
where the sphaleron free energy is found to be 

F SP (T) - ^ (6, 

To one- loop order, the renormalized coupling gn at scale T is expressed in terms of the bare 
(lattice) coupling g as 

Nf is the lattice size in the (Euclidean) time direction. In order to test our lattice approach we 
should verify the quasiclassical prediction for the exponential slope in ^ 

E?±^l~(8LoN t )-^ + const (8) 

if g 2 is small enough. A similar computation for the standard model would be sufficient to 
control at least the order of magnitude of sphaleron transitions close to the electroweak phase 
transition and this would improve our knowledge of this transition rate considerably. 
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In order to determine the sphaleron free energy we shall measure the probability distribution 
of the CS number (after cooling). The CS number itself is given by 

N cs = ^J dxA 1 (9) 

with the vector potential 

A fl = d ll a-sm 2 e/2d^ (10) 

where a is an arbitrary gauge parameter. Convenient choices would be a = for fields close to 
the vacuum at 9 = ir, or a = ip for fields in the upper hemisphere. For a study of the sphaleron 
configurations the latter choice is appropriate. 

In our lattice approach we do not use this expression directly but determine 

A^ = ^(x + d^ X -(d^X + )x) (11) 
through the corresponding CP(l) variables 

A^ is found as the phase of 

X + X2 „ e icosa*/2<W (13) 

where xi an d X2 are on neighbouring sites of link (i. 
The classical equations of motion are given by 

L(x) = d 2 n - lu 2 5 3 = An (14) 

with 

63 = (0,0,1) (15) 

A is a Lagrange multiplier to ensure the constraint. The basic point to observe is that for any 
classical solution L is parallel to n. Therefore, the quantity 

D(x) = L(x) ■ L(x) - [n{x) ■ L(x)} 2 > (16) 

vanishes for any classical solution (including saddle points), and is positive otherwise. The 
integral over D(x) is therefore a convenient new 'action' for defining a cooling procedure which 
does not drive away configurations from saddle points. This procedure can be taken over directly 
to the SU(2) variables (Higgs as well as gauge part) of the electroweak standard model. 

The model was discretized in the standard way on a N t x 128 lattice. This corresponds 
to a temperature 1/N t . We used a heat-bath algorithm to equilibrate the lattice for 36000000 
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sweeps and then a ratio of 8 over-relaxation updates to one heat-bath. We then investigated 10 
blocks of 36000000 sweeps with 20000 configurations measured per block. The averages of each 
of these 10 block measurements were found to be consistent with being statistically independent 
and so were used for an estimate of the statistical error. The mass gap was measured from 
the propagation in the x-direction, namely from the study of the time-slice averages of the 
correlation n\{x)ni{x') + n2(x)ri2(x'). 

To isolate configurations which had a sphaleron, we first required that n^{x) was above 0.5 
on average for a region of x of length 1/lo. For those configurations, we measured Nqs in a 
window from x max — 1/ (2a;) to x max + 1/ (2u>) , where x max corresponded to the maximum value of 
ri3(x). Then configurations with Ncs > 0.45 and n2,{x max ) > 0.90 were classified as sphalerons. 
In order to check the stability of this procedure, we cooled the configurations until the average 
value per link of D was reduced to a fixed value. Then we used the above selection criteria on 
these cooled configurations. The results presented correspond to a reduction of D by a factor 
of approximately 10, although we found that the sphaleron probability was insensitive to this 
threshold value of D (being the same within errors for smaller g 2 ). Our cooling algorithm was 
to replace site variables by an admixture of the original value and that which would minimise 
the action locally. If this procedure reduced D (see above) then it was accepted - otherwise D 
was explicitly minimised which was computationally more demanding since it involved next-to- 
nearest neighbour terms etc. We tried a large number of other choices of sphaleron selection 
criteria - varying thresholds, windows, cooling rate, cooling duration etc and found that the 
changes amounted to an overall constant shift only. Thus the overall normalisation of the values 
for P sp h quoted in the table is not significant, but the dependence on g 2 and N t is. 

Results for the sphaleron probability P sp h per configuration as well as for the mass gap are 
given in Table 1. We used uj = 0.1 in lattice units for Nt = 2 and then varied Nt keeping uoNt 
fixed to check that the results were consistent. Thus it is appropriate to quote the sphaleron 
probability per unit spatial length (in units of uj). Values of P sp h/(128uj) are plotted versus p- 
in Figure 1. Here it is seen that the results for different lattice spacings (ie different uj) are all in 
agreement with each other. Thus we have determined the dependence on g~ 2 of the sphaleron 
production rate per unit length at uj/T = 0.2. This dependence is indeed comparable to the 
quasi-classical expression which is shown by the continuous line (with arbitrary normalisation - 
actually 1.87cT 2 exp(-8uj/(g 2 T)) ). 

The quasi-classical calculation of fluctuations around the classical spaleron gives a result for 
the spaleron probability |25| which is expected to be valid if g 2 T < uj. Our results at g 2 < 0.2 
are thus expected to be approximately reproduced by this approach. As shown in the figure, 
this is indeed the case. 

This close agreement seems a little fortuitous. One reason is that a perturbative expansion in 
terms of lattice parameters does not work very well in the present parameter range, but only at 
much smaller coupling. Thus the measured mass gap and the renormalized coupling should be 
used instead of uj and g 2 in eq.(8). So far, we have not measured the renormalized coupling and 
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Table 1: Monte Carlo results 
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cannot make an estimate of the size of this effect. A somewhat smaller value for the sphaleron 
free energy might also be expected because of the finite range of attraction of the sphaleron in 
our cooling procedure. However, we find very consistent results as the extent of the smoothing 
is varied over a wide range. This suggests that such an effect should be small. 

In conclusion, we have demonstrated a method for measuring Chern-Simons transitions 
which can be carried over to study the bosonic sector of the electroweak standard model. This 
will allow checking various estimates based on resummed perturbation theory for effective po- 
tentials. The outcome is crucial for judging the viability of approaches for generating the baryon 
number of the universe at the electroweak phase transition. 
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Figure 1: The probalitity of a sphaleron-like configuration per unit spatial length as a function 
of 1/g 2 . 
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